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Gravitational waves in cosmology are usually described by transverse traceless tensor perturbations of the back- 
ground Friedmann-Lemaitre- Robertson- Walker (FLRW) metric Q, i.e. 



In a cold matter universe, the linearized gravito-magnetic tensor field satisfies a transverse con- 
dition (vanishing divergence) when it is purely radiative. We show that in the nonlinear theory, 
ii/ , it is no longer possible to maintain the transverse condition, since it leads to a non-terminating 

■ chain of integrability conditions. These conditions are highly restrictive, and are likely to hold only 

in models with special symmetries, such as the known Bianchi and G2 examples. In models with 
realistic inhomogeneity, the gravito-magnetic field is necessarily non-transverse at second and higher 
• order. 
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Q\ • where a 2 jij is the background spatial metric, 7 y Wy = = X> J, Hy, and D is the covariant derivative formed from 
' 7y . An alternative covariant approach due to Hawking is to describe gravitational radiation via the electric and 
q-i magnetic parts of the Weyl tensor, which describes the locally free part of the gravitational field. The advantages 
of this approach are mainly its physical transparency, and the readiness with which one can investigate nonlinear 
j^Tj, extensions of perturbative theory. In particular, using the covariant approach, one can investigate whether the various 
characteristic properties of gravitational waves in the perturbative regime carry through to the nonperturbative level. 
This question is of importance for a covariant understanding of second-order effects in cosmology. 

In general terms, we expect that many linearized features will not extend to second order, precisely because of the 
nonlinearity of the gravitational field, which "creates its own background on which it propagates" . One consequence 
is that the separation of inhomogeneous deviations into scalar, vector and tensor modes is not in general consistent 
beyond the linearized level, and "mode mixing" occurs ||. Indeed, in the fully nonlinear regime, it is probably not 
possible to consistently isolate the radiative part of the gravitational field. Here we investigate whether the transverse 
property can be maintained at second order - and we find that the answer is negative. We show that in general the 
gravito-magnetic field in an irrotational dust universe cannot be transverse in the exact nonlinear theory. This implies 
in particular that at second order, the gravito-magnetic field in cosmology is not in general transverse. Thus there is 
unavoidable "contamination" of tensor perturbations with vector and scalar contributions via the divergence. 



II. COVARIANT APPROACH TO COSMOLOGICAL GRAVITATIONAL RADIATION 



The locally free gravitational field, i.e. that part of the spacetime curvature which is not directly determined locally 
by the energy-momentum tensor, is given by the Weyl tensor C a b c d, which covariantly describes gravitational radiation 
and tidal forces (see e.g. [^],|]-[|). If u a is the cosmic average 4-velocity (with u a u a — —1), then comoving observers 
define the gravito-electric and gravito-magnetic field tensors 

E a b — C ac bdU C U d — E(ab) J Hab = \£acdC cd b e U e = i?( &) j 
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which exemplify a remarkable covariant electromagnetic analogy in general relativity Here e a bc = ilabcdU d is 

the spatial alternating tensor formed from the spacetime alternating tensor rj a b c d, and the angled brackets on indices 
denote the spatially projected, symmetric and tracefree part: 

S(ab) — [h( a c hb) d ~ \h cd h a b\ S c d, 

with h a b = gab + u a Ub the spatial projector into the comoving rest space, and g a b the spacetime metric. 

In linearized theory, the gravito-electromagnetic fields give a covariant description of the gravitational wave back- 
ground in cosmology |^,^|,^||, alternative to the non-covariant description in standard perturbation theory fl], |lo| , |llf . 
A necessary condition for these fields to describe radiation is that their spatial curls should be nonzero; when the 
field is a pure radiation field, additionally their spatial divergences should vanish, in line with analogous properties of 
electromagnetic radiation . The covariant spatial divergence and curl of tensors are defined by [Q 

(divS') a = B b S a b , curl = e cd ( a D c Sb) d , 

where the spatial part D Q of the covariant derivative V Q is given by 

B a A b ... = h a c h b d ---V c A d .... 

Wave-propagation requires non- vanishing curls i.e., 

cmlHab^O^cmlEab, (1) 

if there is gravitational radiation present in the spacetime, allowing the linking of gravito-electric and gravito-magnetic 
variations to each other through the Maxwell-like equations for these quantities, given below. (More precisely, one 
requires curl curl H a b ^ ^ curl curl E a b- In ||, it is shown that a further necessary condition is nonvanishing 
distortions, i.e. D^H^ ^ ^ "DuEt, c \. We will not need these refinements here.) In addition, if the gravitational 
radiation is "pure", i.e., if this radiation encompasses the only deviation of the spacetime (at linear level) from FLRW, 
then the covariant transverse conditions hold: 



(diviT) a = 0= {divE) a . (2) 

This is the simple, transparent and consistent description of perturbative gravitational radiation which flows from the 
covariant approach [|||||. 

In the nonperturbative theory, we can impose the covariant conditions, but wc cannot a priori expect them to 
encode the same simple physical meaning. If we can show that any one of the conditions does not in general hold in 
the nonlinear case, then since these conditions are necessary for pure radiation, that is enough to show a breakdown in 
the linearized description of gravitational radiation at nonperturbative level. Here we focus on the gravito-magnetic 
transverse condition (div H) a = 0, and we show that it is in general inconsistent in the exact nonlinear theory. We do 
not therefore need to investigate the gravito-electric transverse condition (divi?) a = 0, although all the indications 
are that this condition also breaks down at nonperturbative level. 

We consider irrotational dust spacetimes, which provide reasonable models of the late universe, and allow us to 
focus on the purely gravitational dynamics of matter. The known exact solutions with (div H) a = ^ curl H a b are 
spatially homogeneous Bianchi type V, as shown in and G2 solutions with restricted inhomogeneity, found in 
p3[ . In it was argued that in fact (div H) a = does not impose any integrability conditions. Although the basic 
equations and results of fl4|| are correct, there is a flaw in the concluding argument, as we explain below. A simple 
integrability condition implied by (divH) a = was assumed to be satisfied, but it is not identically satisfied in view 
of other valid equations. Rather it imposes non-trivial conditions which have to be checked. Hence, as shown below, 
the transverse condition cannot be obeyed in general, although it can be at the linear level. 



III. INTEGRABILITY OF THE GRAVITO-MAGNETIC TRANSVERSE CONDITION 



The correct statement is that in an irrotational dust universe, (div_ff) a = implies the covariant integrability 
condition 

ee [o-, curl if] = 0, (3) 

where a a b = D(a M b) is the shear, and [A, B] denotes either the tensor commutator of spatial symmetric tensors, i.e. 
[A, B] a b ee A ac B c b — B ac A c b, or its equivalent spatial dual: 
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[A,B] a = e abc A b d B cd = \e abc [A,B] hc . 
To derive Eq. (0) and its consequences, we need some of the covariant propagation and constraint equations, which 



are (12 1 



p = -Qp, (4) 

© = -3 02 -hp-lab^ , (5) 
&ab = — §0CTab — cr c(a (J b) C — -^afc , (6) 

-Eab = -®E ab + 3a c ( a E b ) c + mrlH ab - ^joa-afc , (7) 

H ab = -OH ab + 3a c{a H b) c - curl^ afc , (8) 

(div<7) = |D a 6, (9) 

curl cr ab = i7 afc , (10) 

(div.E) a = [<r,fl] o + ±D p, (11) 

(divlf) = -[cx,£] , (12) 

where p is the energy density, 8 = D a u a is the expansion rate, and a dot denotes the covariant time derivative u a V a - 
(Note that S ab — S/ ab \ since the 4-acceleration vanishes.) 

The general class of irrotational dust spacetimes is characterized by 

p = , u a = = q a , 7T ab = , 

where p is the isotropic pressure, uj a is the vorticity, q a is the energy flux, and it ab is the anisotropic stress. These 
conditions do not lead to any integrability conditions under time evolution, as shown in Q (see also (jl5]-[l8]] ) . By 
contrast, the class of shear-free dust spacetimes, characterized by 

P = , Qa = , (Tab = = TT ab , 

is not in general consistent. Time evolution of the shear-free condition leads to the integrability condition 

e^ a - , 

as shown in [jl9| (see also |2^j2l|]). More severe conditions arise from imposing restrictions on the gravito- 
electric/magnetic fields. 

In linearized theory, (divif) a vanishes by virtue of Eq. (|l2"|), while (div_E) Q = ^D a p, i.e. density inhomogeneity is 
a source for the "Coulomb" part of the gravitational field. In order to isolate the purely radiative part in linearized 
theory, one sets the scalar inhomogeneity D a p to zero. (Note that this does not imply there is no inhomogeneity in 
the matter, only that the inhomogeneity is second order.) Furthermore, in order to remove the remaining scalar and 
vector modes, one sets D a O = 0, which by the constraint equation (^) requires (div a) a — 0. It can then be shown 
that all these conditions corresponding to purely tensor (radiative) modes do not produce any integrability conditions, 
i.e. the covariant description is consistent at the linear level |^,^|. 

At nonperturbative level, (div H) a is no longer identically zero. The vanishing of (div_ff) a is equivalent by Eq. jl^ ) 
to the algebraic condition 

J(°) = [a,E\=0, (13) 

which is of course identically satisfied at the linear level. This covariant commutation property expresses the fact that 
one can choose an orthonormal tetrad which is a simultaneous eigenframe of a ab and E ab . Using the shear propagation 
equation (||) to eliminate E in condition ([l3|), we find that 

[a, &]=Q, 

which allows one to further specialize the eigenframe to be Fermi- Walker transported [Q. Conversely, if the shear 
eigenframe is Fermi-propagated, then Eq. ( |l3| ) necessarily follows, i.e. the gravito-magnetic field is transverse [fl4j| . 

It is possible to investigate consistency of the transverse condition ( |l3| ) using the special tetrad |2^] , but this leads 
to the same conclusion, with no less effort, so we prefer to perform a fully covariant analysis in which the physical 



and geometric quantities themselves are to the fore. The fundamental condition (13) must hold under evolution along 
u a . The first time derivative gives 



3 



X^=-lQX^-^M 0) +^ 1) , (14) 

where X^ is defined in Eq. (||). To derive equation (fL4|), we commute the shear propagation equation (||) with E ab , 
and the gravito-electric propagation equation (Q) with er ab , using basic algebraic properties of the commutator, ft 
follows from equation ( |14| ) that the first time derivative of X^ is not automatically zero by virtue of the original 
condition but yields the integrability condition (|3|). Furthermore, this integrability condition does not follow as 
a consequence of the constraint equations. 

Now the time derivative of the primary integrability condition (JsJ) must also be satisfied. Evaluating this derivative 
requires commuting curl H ab with the shear propagation equation (^), and a a b with (curli? ab )'. There is no basic 
propagation equation for curliif ab . However, the propagation of curlif ab is indirectly determined by the curl of the 
gravito-magnetic propagation equation (|^), together with the identity ]T^ | 

(curlS ab )' = curlS ab - ±6curlS ab - <r e c e cd{a D e S b) d + 3H c{a S b) c . (15) 

Then using Eqs. (H) and (15), we find that 

iM = -20iP-a a b lU-lW, (16) 

where 

X a 2 ) = [a, curl curl E] a - 3[a, curl {aH)] a + [E,cuilH] a + 3H a b [a,H] b 
+ <j bc H bc D a e - (H ac <j cb + \a ac U cb ) D b 6 - a b d <r bc D c H ad 

+ {(Ta C VM + \°ah<J C d) V C H M . (17) 

Here (AB) ab = (AB)^ = A c ( a B b ^ c , and we used the properties e a bc£ de ^ = Slh^h^h^ and e a bc£ dec = 2h[ a d h b f. 
We can rewrite this expression for by using the identity [h2| 

curl curl S ab = -D c D c ^ Qb + |^ (a (div5) b> + (p- ±0 2 ) S ab 

+ 3S c {a {Eb) C — |0 f7 &> C } + &cdS C (Tab — S cd (T ca (J bd + <J cd <J c ( a S b )d ; 

and the identity |^| 

D c S a b = D^Saty — |e dc ( a curlS' b ) d + |(div S)^ a h b ^ c . 

However the expression remains intractable. 

Clearly satisfying the primary integrability condition (|3|) on an initial surface does not in itself guarantee that 
remains zero. Instead, the vanishing of requires, by Eq. (|l6|), an independent secondary integrability condition 

X (2) =0, (18) 

which is of course identically satisfied at the linear level. The form of Eq. (|l7|) shows that the integrability conditions 
grow more complicated at each stage, involving a growing order of spatial derivative as well as increasing algebraic 
complexity. The constraint equations do not lead to any simplification of this integrability condition. 

It is clear that further time evolution will produce a third, and then a fourth, fifth, . . . integrability condition. 
Each condition is more complicated than its predecessor, and is not identically satisfied in general by virtue of earlier 
conditions or the constraint equations. Furthermore, eliminating terms higher than second order does not change this 
fundamental feature, although it does simplify the integrability conditions somewhat. At second covariant order, the 
primary integrability condition (||) is unchanged in form, but the secondary integrability condition ( |l8| ) becomes 

X {2) = [a, curl curl E] + [E, curl H] = . (19) 

The next integrability condition will involve third order spatial derivatives, and so on. 

IV. CONCLUSION 

We have shown that the gravito-magnetic transverse condition (diviJ) a = 0, taken over from linearized theory 
into the nonlinear theory, leads to a chain of derived integrability conditions, the first two of which are Eqs. (0) and 
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(|l8|). This result corrects the assertion in |Q that condition (^) (and thus its derivatives) is identically satisfied. The 
logical flaw in the argument presented in pTf is essentially the implicit assumption that the propagation equations 
@ _ ® are identically satisfied. We have corrected that error here, by imposing the propagation equations as effective 
conditions on the evolution of (divH) a — 0, equivalently [cr,E] = 0. 

The consequence of this is that in general, the gravito -magnetic field cannot he transverse at second (and higher) 
order, so that it cannot be purely radiative beyond linear order. In the fully nonlinear case, when one cannot really 
isolate gravitational radiation, our result implies that the transverse condition is in general inconsistent, i.e. generic 
inhomogeneity within the class of irrotational dust spacetimes is inconsistent with a transverse gravito-magnetic field. 
These results are independent of gauge choices or of coordinates, since they are fully covariant. 

There are however special cases when the integrability conditions are satisfied, the most important of which is 
the case of almost-FLRW universes. The condition (^) and all its derivatives are satisfied to first order in this case, 
which simply reflects the fact that the covariant formulation of the transverse condition for perturbative gravitational 
radiation is consistent at the linear level. 

The integrability conditions are also satisfied if H a \, = 0, which defines the "silent" models. However, in that case, 
there is another chain of integrability conditions flowing from the vanishing of H a b at all derivative levels. As shown in 
^,^4|, these integrability conditions are even more severe than in the (div_ff) a = case. The H ab = models are thus 
in general inconsistent (see also [p5| ), and unlikely to extend beyond the known special cases where the integrability 
conditions collapse to identities. These special cases are some Bianchi and Szekeres solutions (see |2^,Q). Such a 
chain of derived integrability conditions of growing complexity exists also in the "anti-Newtonian" models E ab = 
(see p7|), where the only known (irrotational dust) solution is the degenerate FLRW solution. It is clearly very 
restrictive to impose conditions on the gravito-magnetic and -electric fields of irrotational dust spacetimes. 

In the (divH) a = =^ H ab case, among the known special exact solutions are Bianchi type V jl2| and diagonal (in 
comoving coordinates) G 2 solutions jl3). In these solutions, the integrability conditions are identically satisfied by 
virtue of the special structure of curl_ff a h. If we choose a Fermi-propagated shear eigenframe (as described above), 
then in these solutions, car\H a b remains diagonal at all times, so that the primary integrability condition (||) is 
identically true, as are all its derivatives. For example, for G 2 solutions of the form 

ds 2 = -dt 2 + A 2 (t,x)dx 2 + B 2 (t,x)dy 2 + C 2 {t,x)dz 2 , 

in coordinates that are comoving with the fluid 4- velocity u a , it follows that 

which leads to a diagonal curl_ff a &. The Bianchi type V solution has cuvlH a i, cx a a \, in the shear eigenframe. In 
the Bianchi case, it is possible to identify further examples pti|] . The type V case is distinguished as the only 
(div H) a — 2^ Hab solution in Bianchi class B, and it has (div E) a ^ 0, so that it is not purely radiative in the sense 
of condition (g). (Note that the diagonal G 2 solutions also have (divE) a 7^ 0.) 

Indeed it is clear from our analysis that a similar chain of integrability conditions will flow from the gravito-electric 
transverse condition in Eq. (||), which is equivalent to 

lB a p+[a,H] a = 0. 

The independent chain of conditions arising from repeated time differentiation of this condition reinforces the unlike- 
lihood of finding inhomogeneous transverse radiative solutions beyond the linear level. 

The Bianchi class A solutions all have vanishing divergence for both E ab and H ab , and so provide spatially homo- 
geneous pure radiative examples, except for types I and LRS VIIo, which are in the degenerate silent case H ao = 0. 
Thus in Bianchi class A spacetimes it is possible to find exact solutions satisfying both the radiative transverse con- 
ditions ^) and the propagating radiative conditions We give here two such examples, using the formalism and 
notation of p0| . In a shear eigenframe {e Q } (where a = 1, 2, 3), a tracefree diagonal tensor S a p has curl 

curlS* Q/ 3 = e 1& ( a d s S l3 ) 1 + \S a pn 5 s - 3n 7 (ct S , (3 ) 7 + 8 a pS lS n< s , (20) 

where d a is the directional derivative along e Q , and the tetrad commutation quantities may be diagonalized: n a p = 
diag(n l7 n 2 , n 3 ). Using the irreducible components S + ~ — §Sii and 5_ = ^v / 3(5 , 22 — £33), the Bianchi type II 
solutions (characterized by n\ > and n 2 = = 713) have nonvanishing components of cur\S a p 

curl 5+ = — |ni5+ , curl SL = ^niS- , 

and nonvanishing components of curl curl 
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curl curl 5+ = |(ni) a 5 + , curlcurlS_ = ^(ni) 2 ^ . 
The type VIo solutions with n a a = (and n-i — —n\ < 0, 77,3 = 0), have nonvanishing components of cuilS a p 

curlS+ = -±V3m(\/3S + - S-) , curlSL = ±V3m(S+ + VSS-) , 
and nonvanishing components of curl curl S a j3 

curl curl S± — 3(ni) 2 S± . 

The tensor S a p in both examples may then be chosen to be either the gravito-electric tensor E a p or gravito-magnetic 
tensor H a p. 

The importance of these examples is twofold: (a) they show that exact purely radiative solutions can exist in 
the minimal covariant sense, and (b) they show that the Petrov type N characterization of a radiative field is more 
restrictive, because there are no Petrov type N dust metrics in the full nonlinear theory [^8| , ^9| . However, because 
these examples are spatially homogeneous models, information is not being conveyed from place to place by the waves; 
these are "standing" waves rather than traveling (propagating) waves. By contrast, in the linearized case, arbitrary 
information can be conveyed by the gravitational waves || . 

It is possible that further, inhomogeneous exact solutions can be found, but they would also be very special cases, 
in order to satisfy the chain of integrability conditions. Realistic inhomogeneous models cannot be expected to satisfy 
these conditions. 
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